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Abstract 

The methane-nitrogen phase diagram of Prokhvatilov and Yantsevich (1983) 
indicates that at temperatures relevant to the surfaces of icy dwarf planets 
like Pluto, two phases contribute to the methane absorptions: nitrogen satu¬ 
rated with methane N 2 :CH 4 and methane saturated with nitrogen CH 4 :N 2 . 
No optical constants are available so far for the latter component limiting 
construction of a proper model, in compliance with thermodynamic equilib¬ 
rium considerations. New optical constants for solid solutions of methane 
diluted in nitrogen (N 2 :CH 4 ) and nitrogen diluted in methane (CH 4 :N 2 ) are 
presented at temperatures between 40 and 90 K, in the wavelength range 
1.1-2.7 /im at different mixing ratios. These optical constants are derived 
from transmission measurements of crystals grown from the liquid phase in 
closed cells. A systematic study of the changes of methane and nitrogen solid 
mixtures spectral behavior with mixing ratio and temperature is presented. 
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1. Introduction 


Pluto, Eris, and Makemake, unlike most trans-Neptunian objects (TNOs) 
with water-ice rich or featureless surfaces (Barucci et ah, 2008), display in¬ 
frared spectra dominated by methane ice (Brown, 2008). These three TNOs 
are often compared with Neptune’s large satellite Triton, since its spectrum is 
dominated by methane ice and it is thought to have formed similarly to Pluto, 
Eris, and Makemake, prior to its capture into a retrograde orbit around Nep¬ 
tune. In addition to methane ice, nitrogen ice has been directly detected on 
Pluto and Triton via the 2.148-pm absorption band (Cruikshank et ah, 1984; 
Owen et ah, 1993). Additional evidence for nitrogen ice on Pluto and Triton 
comes from shifts of their methane absorption bands to shorter wavelengths, 
which occurs, as described by Quirico and Schmitt (1997), when methane 
is dissolved at low concentrations in a matrix of solid nitrogen. Moderately 
high signal-to-noise spectra of Makemake and Eris show no evidence for the 
presence of the 2.148-/im nitrogen absorption feature. Their methane bands 
do present subtle shifts, albeit smaller than the shifts measured in spectra 
of Pluto and Triton. A laboratory study by Brunetto et ah (2008) showed 
that smaller shifts correspond to higher methane abundances. The lack of 
the 2.15-/im nitrogen absorption band and the smaller methane wavelength 
shifts led several authors (Brown et ah, 2007; Alvarez-Candal et ah, 2011; 
Merlin et ah, 2009) to the conclusion that Eris and Makemake are not nitro¬ 
gen dominated, contrary to Pluto and Triton. 

Thermodynamic equilibrium dictates that if methane and nitrogen ices 
are both present, for most of the range of possible nitrogen/methane relative 
abundances, two distinct phases must coexist at temperatures relevant to the 
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surfaces of these icy dwarf planets (Prokhvatilov and Yantsevich, 1983; Lu- 
nine and Stevenson, 1985): methane ice saturated with nitrogen and nitrogen 
ice saturated with methane. 

Tegler et al. (2010) demonstrated that the depth of the methane and nitro¬ 
gen absorption bands and the wavelength shift of the observed methane ab¬ 
sorption features should not be used as proxy for the methane-nitrogen mix¬ 
ing ratio. The phase composition as dictated by thermodynamic equilibrium 
must be taken into account. In particular, Tegler et al. (2010) model each of 
the observed methane absorption bands with a binary mixture of methane 
ice saturated with nitrogen and nitrogen ice saturated with methane, accord¬ 
ing to the methane-nitrogen phase diagram of Prokhvatilov and Yantsevich 
(1983). This technique, applied so far to the cases of Eris and Pluto only, 
results in the Ending that the bulk volatile composition of Eris is similar to 
that of Pluto, with both objects being dominated by nitrogen ice. A correct 
modeling of the methane absorption bands has strong implications not only 
on the methane-nitrogen mixing ratio, but also enables exploration of stratifi¬ 
cation properties as well as heterogeneity of these targets (Tegler et ah, 2010, 
2012). Because of the lack of absorption coefficients for methane saturated in 
nitrogen and for nitrogen saturated in methane at the appropriate tempera¬ 
tures in the visible wavelength ranges, Tegler et al. (2010) approximate the 
methane-dominated phase by using pure methane absorption coefficients and 
the highly diluted phase by shifting pure methane coefficients by amounts 
seen for highly diluted samples (Quirico and Schmitt, 1997). 

In this paper we provide optical constants in the wavelength range 1.1- 
2.7 pm of solid solutions of methane diluted in nitrogen, N 2 :CH 4 , and ni- 
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trogen diluted in methane, CH 4 :N 2 , at temperatures between 40 and 90 K 
and at different mixing ratios (http://www 2 .lowell.edu/users/grundy/ 
abstracts/2015.CH4+N2.html), allowing a proper model to be constructed 
for any TNO where the methane/nitrogen ratio falls between the two sol¬ 
ubility limits such that both saturated phases are present. Gaining more 
detailed knowledge of the methane/nitrogen mixing ratio and phase state 
of Pluto, Triton, Eris, and Makemake will enable a better understanding of 
the processes responsible for volatile loss and retention on TNOs (Schaller 
and Brown, 2007). It will also help constrain the seasonal behaviors of their 
atmospheres, supported by vapor pressure equilibrium with surface ices. It 
will also shed light on the photolytic and radiolytic chemistry that can occur 
within the surface ices, since radicals produced by energetic radiation will 
encounter and react with different molecules within a nitrogen-dominated or 
a methane-dominated solid phase. 

2. Laboratory Experiments 

The experiments reported here were conducted in a new laboratory ice 
facility located in the Department of Physics and Astronomy of Northern 
Arizona University. A detailed description of this facility is given by Tegler 
et al. (2010, 2012 ) and Grundy et al. (2011). We used the closed cell technique 
for the ice sample preparation. This technique consists in the growth of 
crystals from the liquid phase in a closed cryogenic cell. Grystals are grown 
as follows. The sample is prepared in gas form in a 2 liter mixing volume. 
Here we report experiments obtained mixing methane (GH 4 ) and nitrogen 
(N 2 ). The purities of the gases used are 99.999% for GH 4 and 99.9% for N 2 , 
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as reported by the vendors. The gas is set to flow into an empty cell, which 
is at a temperature slightly higher than the melting point of the ice sample. 
In the case of mixtures, the melting point depends on the mixing ratio of the 
gases and we use, as reference, the CH 4 -N 2 phase diagram of Prokhvatilov 
and Yantsevich (1983). The gas, once in the cell, condenses immediately to 
liquid. The liquid is frozen by reducing the temperature in the cell at a rate 
of 0.1 K minute"^. A thermal gradient is maintained within the cell with the 
top and bottom heaters such that the crystal grows from the bottom to the 
top. Once the sample is frozen, the vertical thermal gradient is minimized. 
After the initial ice spectrum is recorded, the temperature is ramped down 
at 0.1 K minute“^. We recorded spectra at temperatures between 40 and 
90 K. 

A 5 mm cell, with sapphire (AI 2 O 3 ) windows, was used (Figure 1). Thin¬ 
ner samples, needed in the case of mixtures with high CH4 content, were 
prepared by use of a transparent potassium bromide (KBr) or calcium flu¬ 
oride (CaF 2 ) spacer between the windows. These KBr and CaF 2 spacers 
occupied only part of the cell, enabling the choice of the full or reduced 
thickness, simply by moving the sample relative to the spectrometer beam. 
The thickness of the sample in absence of spacer is equal to the cell depth, 
d = d', which is known. On the other side of the cell, due to the uncertain¬ 
ties on the depth of the KBr and CaF 2 spacers, the sample has an unknown 
thickness, d = d". This is computed with a liquid CH4 experiment. Two 
transmission spectra are recorded, T' and T", by positioning the spectrome¬ 
ter beam away and in correspondence of the spacer, respectively. From the 
two transmission spectra, we compute the imaginary part of the refractive 


5 


/,• 


mo 






t 

mce («J0) 


! 

d\ ^ 

i 


d 


mKBr 


d\+d2=d 

\ 

mice (mo) 


i 



mAho, 




mo 


filled i^empiy^ 


Figure 1: Schematic view of the cell interior. See text for details. The incident radiation li 
is transmitted through media of different refractive index m. lempty and Ifuied represent 
the transmitted radiation when the cell is empty (the vacuum toq is between the sapphire 
windows iriAhOa, and the potassium bromide spacer rriKBr) and filled by the sample 
{mice), respectively. 

index, k' and k”, which will be function of d' and the unknown d", respec¬ 
tively. Because the imaginary part of the refractive index, k, is a property 
intrinsic to the material, we compute d" by setting k'=k'' (see Section 3). 

Spectra were recorded with a Nicolet Nexus 670 Fourier transform in¬ 
frared (FTIR) spectrometer, covering the range 2799.9 - 12000.1 cm“^ at a 
sampling interval of 0.24 cm“^, resulting in a spectral resolution of 0.6 cm“^ 
(FWHM of unresolved lines). We averaged over 100 spectral scans to improve 
the signal-to-noise ratio. 
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Figure 2: The transmission spectra as a function of wavelength and wavenumber, recorded 
with the cell empty (dashed line) and filled (continuum line) by pure CH 4 ice at 40K. The 
upper spectrum has been shifted along the y-axis by 0.4 for clarity. 


3. Data Analysis 

3.1. Calculation of n and k 

The purpose of our measurements is to compute optical constants of 
methane-nitrogen ice mixtures, CH 4 -N 2 . The optical constants of a mate¬ 
rial are the real and imaginary part of the complex refractive index m{u) = 
n(z/) -b and are both functions of the frequency u. Let T{v) be the 

transmission spectrum of an ice sample within a cell of thickness d (Fig¬ 
ure 1 ), we compute the imaginary part of the refractive index, k, via the 
Beer-Lambert absorption law: 


k{u) = 


Aiivd 


lnT(z/) 


( 1 ) 


In order to remove the flux distribution of the illumination source, the trans¬ 
mission function of the spectrometer and the cell, the detector spectral sen¬ 
sitivity function, and the water vapor contamination (Grundy and Schmitt, 
1998; Tegler et al., 2012), we compute the transmission spectrum T(z/) as: 

Tfillediy^ 


T{u) = 


Tempty {y^ 


( 2 ) 
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Figure 3: Top, middle, and bottom rows show the transmission spectrum, the imaginary 
part of the refractive index, and the real part of the refractive index of pure CH 4 at 40K, 
respectively, in four different bands. 

where Tfmediy) and Tf^mptyW) represent the transmission through the cell 
when filled by the sample and empty, respectively (Figure 2). Specifically, 
Tempty{i^) IS the average of the transmission spectra through the empty cold 
cell recorded before and after the ice sample experiment. The transmission 
spectra, whether the cell is filled or empty, present a high frequency interfer¬ 
ence pattern due to the windows in the cell, which we remove using a Fourier 
filter. The resultant transmission spectrum T(z/) is affected by subtle slopes 
from a variety of sources (Grundy and Schmitt, 1998). In order to remove 
these slopes, we divide the spectrum by the fit to the continuum regions ad¬ 
jacent to the absorption bands to be analyzed. The continuum of our spectra 
is well fit by a third order polynomial. 

Given the imaginary part of the refractive index, k{u), we calculate n{u) 














































from the Kramers-Kronig dispersion relation 


n(z/) = no + - 

TT 


u'k{u') 


z/' 


/2 


-dv' 


( 3 ) 


where uq is the value of n at the high-frequency end of the infrared, and the 
integration is over the infrared region of the spectrum. In the case of pure 
CH 4 , uq is equal to 1.32 at 15800 cm“^ (Trotta, 1996). We will adopt the 
same value also for CH 4 -N 2 mixtures. 

After computing n(z/) and we rehne these calculations by taking into 
account hrst-order differences in transmission through boundaries between 
cell windows and vacuum in the case of empty cell, Temptyiy)) and between 
cell windows and ice in the case of hlled cell, Tjuiedii^)- If the spectra are 
acquired by putting the beam on the side of the cell containing the KBr or 
CaF 2 , we have to take into account the interfaces between the ice or vacuum 
and the KBr or CaF 2 , as well (Figure 1). In this case, we have: 


k{v) = 


TAl203,ice Tice,KBr TKBr,ice Tic£,Al203 lnT'(z/) 


dvrz/d L Ta 12O3,0 To,KBr TKBr,0 Tq,a1203 


( 4 ) 


where Tpq represents the transmission across the boundary between the media 
with refractive index rup {irtp = Up + ikp) and rUg. Multiple reflections that 
are set up within the slab as the wave bounces back and forth at the top and 
bottom interfaces are considered negligible. For the case of normal incidence 
we have: 


T = 

-^pq 


4 n, 


■pq 


(jlpg + 1)2 kl, 


pq 


where Upg and kpg are given by: 


( 5 ) 
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+ kl 


( 6 ) 
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The transmission spectrum and optical constants (n and k) of pure CH 4 ice 
at 40 K are shown in four different bands in Figure 3. 


3.2. Pure CH 4 ice 

Grundy et ah ( 2002 ) present an exhaustive analysis of the temperature- 
dependent near-infrared absorption spectra of pure CH 4 ice, available at tem¬ 
peratures between 30 and 93 K. Our measurements for pure CH 4 ice at 40 K 
are compared with that by Grundy et ah (2002) in Figure 4. The absorp¬ 
tion spectra, a, by Grundy et ah ( 2002 ) have been converted for comparison 
purposes into the imaginary part of the refractive index, fc, by means of the 
dispersion relation a = Airk/X. While the data by Grundy et ah ( 2002 ) are 
obtained combining measurements acquired with different cell thicknesses, 
enabling measurement of the weakest and strongest GH 4 absorptions, our 
data for pure GH 4 ice correspond to a single sample of thickness 0.0252 cm. 
This prevents us from measuring absorption where it approaches zero or 
100%. Therefore, we provide our results in narrow blocks of wavelengths 
covering regions of intermediate absorption. Differences in band strength 
between our measurements and those by Grundy et ah ( 2002 ) are observed, 
as shown in Figure 4 for T= 40 K. Similar discrepancies between the two 
data sets are observed at all temperatures. This mismatch could be due to 
continuum removal, uncertainty in the sample thickness and/or temperature, 
transmission across the boundaries not taken into account when converting 
the absorption coefficients reported by Grundy et al. ( 2002 ) into k. 
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Figure 4: Comparison between the imaginary part of the refractive index of pure CH 4 ice 
presented in this work (red solid line) and by Grundy et al. (2002) (black dash dot line) 
at 40 K. The ratio between the two data sets is shown in the bottom panel. 

The general character of the temperature dependence of the pure CH 4 
ice absorption coefficients discussed by Grundy et al. ( 2002 ) is confirmed 
by our data (Figure 5): lower temperature spectra exhibit narrower and 
higher peak absorptions which broaden and decrease in strength as the ice 
is warmed. Most of the CH 4 ice bands do not show significant temperature- 
dependent wavelength shifts. Nevertheless subtle shifts occasionally arise; a 
good example can be seen at 6735 cm“^. Figure 6 displays the peak position 
of the z/i-|- 3 i /4 CH 4 ice band centered around 6735 cm“^ as a function of 
temperature. The peak position has been estimated as the average of the 
position where the maximum of the band occurs and the maximum of the 
gaussian function fit to the data around the band center. The error is the 
dispersion between the two estimates. 
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Figure 5: k spectra of pure CH 4 ice are shown for temperatures between 33 and 90 K. 
The inset panel is an enlargement of the region around 5800 cm“^. In order to show, in 
greater detail, the thermal broadening of the k spectra of pure CH 4 ice with temperature, 
the latter have been normalized to their maximum value in the range from 5750 cm“^ to 
5850 cm-i. 
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Figure 6 : Panel (a) shows k spectra of pure CH 4 ice at temperatures between 33 and 90 K 
in the range between 6700 and 6810 cm“^. The spectra have been normalized to their 
maximum value. The peak position as a function of temperature for the izi+ 3 iz 4 CH 4 ice 
band is shown in panel (b). 
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3.3. CH 4 -N 2 Solid Mixtures 


The methane and nitrogen stoichiometry of a CH4-N2 solid mixture is 
estimated by dividing the integrated k spectrum of the mixture by the corre¬ 
sponding integrated k spectrum of pure CH4 ice. This results in the fraction 
of CH4 in the mixture. Given the temperature dependence of the pure CH4 
ice optical constants (see Section 3.2), k spectra of CH4-N2 mixtures and pure 
CH4 acquired at the same temperature are compared. This procedure is ap¬ 
plied in the frequency ranges from 5300 - 6440 cm“^, 6700 - 7650 cm“^, and 
8100 - 9000 cm“^. Fcha,^ in Table 1 represents the CH4 fraction computed in 
the interval centered around v. The values obtained at all temperatures and 
frequency ranges for a given CH4-N2 mixture are averaged, resulting in the 
CH4 stoichiometry of our sample {CH 4 [%]). The mean CH4 content, CH 4 
[%], is listed in Table 1 for all our experiments. The error on the CH 4 con¬ 
centration is given by the scatter of the values computed at all temperatures 
and frequency ranges and it is generally less than 5%. The measurements we 
present encompass both sides of the CH4-N2 phase diagram of Prokhvatilov 
and Yantsevich (1983, Figure 7) enabling us to characterize the spectral be¬ 
havior of solid solutions of CH4 diluted in N2 (N2:CH4) and N2 diluted in 
CH4 (CH4:N2) at different temperatures and mixing ratios. Example opti¬ 
cal constants of CH4:N2 (92.04% CH4) and N2:CH4 (12.48% CH4) at two 
different temperatures at low spectral resolution are reported in Table 2. 

Laboratory data of CH4:N2 and N2:CH4 at the solubility limits of CH4 and 
N 2 in each other, which we indicate with CH 4 :N 2 and N 2 :CH 4 respectively, 
give the means to generate optical constants of CH4-N2 mixtures with any 
bulk composition between the solubility limits, at a given temperature (Fig- 
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Figure 7: The CH 4 -N 2 binary phase diagram generated from X-ray diffraction studies 
by Prokhvatilov and Yantsevich (1983) is shown. The symbols indicate the different set 
of transmission measurements performed to conduct a systematic study of the changes 
in CH 4 -N 2 solid mixtures spectral behavior with mixing ratio and temperature. Mea¬ 
surements of CH 4 :N 2 and N 2 :CH 4 solid solutions acquired at different temperatures are 
represented by the same color and symbol. Between the solubility limits we have a mixture 
of two solid solution phases: CH 4 ice saturated with N 2 , CH 4 :N 2 , and N 2 ice saturated 
with CH 4 , N 2 :CH 4 . 


ure 7). This is possible through the lever rule, which is a tool used to deter¬ 
mine the molar percentage of each phase (in our case CH 4 :N 2 and N 2 :CH 4 ) in 
a two phase alloy in equilibrium. The molar percentage of the CH 4 :N 2 phase 


in a CH4-N2 mixture is 




Xch4 


qCHa 

‘^N7:CH4 


0CH4 

‘^N7:CH4 


0CH4 

‘^CH4;;N2 


( 8 ) 


where and represent the solubility limit of CH4 in the 
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CH 4 :N 2 and N 2 :CH 4 phases, respectively, and Xch^ is the overall CH 4 abnn- 
dance in the mixtnre. At 55 K the solnbility limits of CH 4 in the CH 4 -rich 
and N 2 -rich satnrated phases are close to 82.69% and 16.35%, corresponding 
to measnrements of optical constants we made in the laboratory (hlled dia¬ 
mond and open downward triangle in Fignre 7). It is therefore possible to 
generate synthetically the optical constants of any CH 4 -N 2 mixtnre at 55 K. 
To a hrst order approximation onr data behave consistent with expectation 
from the lever rnle (see Fignre 8 ). 


^ Wavelength (gm) 

^ 1.82 1.79 1.75 1.72 1.69 1.67 1.64 1.61 1.59 



Figure 8 : k spectrum of a CH 4 -N 2 mixture measured in the laboratory at 55 K and 
containing 71.88% of CH 4 (black solid line) and the synthetic fc-spectrum (dashed red 
line) generated applying the lever rule to the optical constants of CH 4 :N 2 (dash-dot 
orange line) and N 2 :CH 4 (dash-dot blue line) phases at 55 K. The synthetic spectrum 
reproduces the measured one well. 

Methane-nitrogen mixtnres present a temperatnre-dependent behavior 
similar to that of pnre CH 4 (Section 3.2). We characterize the thermal 
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Figure 9: Full width at half maximum (FWHM) as a function of temperature and stoi¬ 
chiometry for three isolated CFl 4 -ice bands. Samples with different CH 4 and N 2 content 
are represented by different symbols (see legend in Figure 7). The dashed and solid lines 
represent the best fit of the CH 4 :N 2 and N 2 :CH 4 data, respectively. 

broadening of CH 4 ice absorption bands by measuring their full widths at 
half maximum (FWHM) as a function of temperature and stoichiometry 
(Figure 9). This analysis has been conducted for those CH 4 -bands that are 
isolated, not influenced by the wings of adjacent bands. These are the 1 ^ 2 + ^' 3 , 
^^2 + t '3 + i^ 4 , and i/ 2 + 2 z /3 CH 4 -band, centered, in a sample of pure CH 4 at 90 K, 
at 4526, 5800, and 7487 cm“^, respectively. All bands we measured show a 
nearly linear dependence of FWHM with temperature, at all stoichiometries. 
Furthermore, the CH 4 bands of CH 4 :N 2 samples broaden faster than those 
of N 2 :CH 4 , as expected given the higher CH 4 content in CH 4 :N 2 than in 
N 2 :CH 4 mixtures. 

The imaginary part of the refractive index, fc, of CH4-N2 ice mixtures at 
different mixing ratios at T=60 K is shown in Figure 10. CH4 when dis¬ 
solved in N 2 presents absorption bands shifted toward shorter wavelengths 
compared to the central wavelengths of pure CH4. This shift varies with 
the CH4 abundance in the mixture: the larger the CH4 concentration the 
smaller the blueshift. We estimate the spectral shift by means of a cross 
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correlation technique. We shift a band in the pure CH4 spectrum relative to 
the same band in the mixture spectrum and for each shift a is computed. 
The versus shift behavior is htted with a parabola, whose minimum po¬ 
sition corresponds to the shift solution for the analyzed band. This method 
is described in detail by Tegler et ah ( 2008 ). Another possible approach is, 
for a given band, to determine the difference between the peak positions in 
the CH4-N2 and pure CH4 samples, with the peak position estimated as de¬ 
scribed in Section 3 . 2 . The solutions obtained with these different techniques 
are averaged, and the dispersion between them is taken to be representative 
of the uncertainty. In order to correctly determine the blueshift of the CH4- 
N2 bands relative to pure CH4, spectra acquired at similar temperatures are 
compared, given the temperature dependent behavior of the CH4 bands peak 
position (Section 3 . 2 ). The blueshift of four CH4 bands versus CH4 abun¬ 
dance is reported in Figure 11 . For a given CH4 stoichiometry, different mea¬ 
surements correspond to different temperatures. No signihcant variations in 
the wavelength shift with respect to temperature are observed. Filled right¬ 
facing triangles in Figure 11 , contrary to the other symbols that represent 
samples generated in the laboratory (see legend in Figure 7 and Table 1 ), 
refer to optical constants computed synthetically by applying the lever rule 
at 55 K for a range of CH4 abundance between 30 % and 70 %, in 10 % steps. 
The N2-enriched and CH4-enriched sides of the phase diagram present a dif¬ 
ferent linear trend of the blueshift versus CH4 abundance (dashed and solid 
lines in Figure 7 ). Such trend varies from band to band. 
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Figure 10: Top panel: The imaginary part of the refractive index, k, of CH 4 -N 2 ice 
mixtures at different mixing ratios at T = 60 K over the range from 5490 to 6300 cm“^. 
Bottom Panel: Expanded view of the band. 

4. Applications to Pluto 

Efforts to quantitatively model Pluto’s spectrum over the course of the 
years identihed several important issues. The hrst concerns the state of 
CH 4 ice on Pluto’s surface. Seen at high spectral resolution, the CH 4 bands 
in Pluto’s spectrum are shifted toward shorter wavelengths compared to 
the central wavelengths of pure CH 4 obtained in the laboratory, implying 
CH 4 being dissolved at low concentrations in a matrix of solid N 2 (Schmitt 
and Quirico, 1992; Quirico and Schmitt, 1997). Additionally, the identih- 
cation of the 1.69-/im band in Pluto’s spectrum has been interpreted so far 
as the evidence of the presence of pure CH 4 on Pluto. Indeed, the 1.69- 
/im CH 4 feature has never been observed in any sample of CH 4 diluted at 
low concentrations in N 2 (neither in a- or /3-N2 phases) at any temperature. 
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Figure 11: Blueshifts of the V 2 + Vd, (4526 cm“^), V 2 +V 3 + va (5800 cm“^), 21/3 (5991 cm“^), 
and V 2 + 21/3 (7487 cm“^) bands versus CH 4 abundance. Symbols and colors correspond 
to samples with different CH 4 and N 2 stoichiometries as described by the legend in Fig¬ 
ure 7 (see text for details). Solid and dashed lines are the linear fit to the shift versus 
CH 4 abundance for percentages greater and less than 50%, respectively. 

These considerations together with the lack of optical constants for CH 4 with 
a small fraction of N 2 are the main argnments behind the approach taken 
in analyzing Pinto spectra in the last decade {e.g., Doute et ah, 1999; Olkin 
et ah, 2007; Protopapa et ah, 2008), employing pnre CH 4 and CH 4 diluted 
at low concentrations in N 2 . 

The set of optical constants presented in this work enable us to compare 
and contrast samples with CH4 dissolved at low concentrations in a matrix of 
solid N2 (N2:CH4), and CH4 with a small fraction of N2 (CH4:N2). These data 
reveal that the 1.69-/im band of Pluto’s spectrum can no longer be considered 
as evidence for the presence of pure CH4 on Pluto’s surface, as this feature 
is observed not only in samples of pure CH4 (at temperatures below ~60 K, 
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Figure 5), but also in CH 4 :N 2 samples (Figure 10). Also, the wavelength shift 
of Pluto’s CH 4 bands with respect to that of pure CH 4 , does not necessarily 
imply the presence of CH 4 being dissolved at low concentrations in a matrix 
of solid N 2 , as the wavelength shift is a linear function of the CH 4 abundance 
(Figure 11). 

The quantitative analysis of the wavelength shift versus abundance indi¬ 
cates that the wavelength shift of the CH4 bands is indeed a good dilution 
indicator. The linear hts of the wavelength shift as a function of CH4 abun¬ 
dance for different CH4 bands provide the tools to infer the CH4-N2 mixing 
ratio on the surface of Pluto, or any TNO covered by CH4 and N2. Also, it is 
extremely important to stress that the linear dependency of the shift versus 
CH4 abundance changes from band to band. This needs to be taken into 
account when comparing the shift of deeper and shallower bands to address 
stratihcation properties (Merlin et ah, 2009). 

Our new optical constants enable us to model for the first time Pluto’s 
near-infrared spectroscopic measurements in compliance with thermodynamic 
equilibrium, superseding more than a decade of work in which Pluto’s spec¬ 
trum was interpreted in terms of a combination of “diluted” and “pure” 
CH 4 components. To account for thermodynamic equilibrium on Pluto’s sur¬ 
face, we need instead to consider methane saturated with nitrogen, CH 4 :N 2 , 
and nitrogen saturated with methane, N 2 :CH 4 , as dictated by the N 2 -CH 4 bi¬ 
nary phase diagram (Figure 7, Prokhvatilov and Yantsevich, 1983) and sug¬ 
gested by Trafton (2015) and Cruikshank et al. (2015). Specihcally, for a 
plausible Pluto volatile ice temperature of 40 K (Tryka et ah, 1994), the dia¬ 
gram shows that the solubility limits of CH 4 and N 2 in each other are about 
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5% (N2:CH4, with 5% CH4) and 3% (CH 4 :N 2 with 3% N 2 ). While we have 
optical constants available for the former phase (see Figure 7), the latter is 
not available yet. However, we have laboratory data covering a wide range 
of temperatures and mixing ratios to synthetically reproduce the lacking set 
of optical constants. The laboratory measurements presented in this work 
cover the wavelength range of the LEISA spectrometer on board of the New 
Horizons spacecraft (Young et ah, 2008), which will flyby Pluto in July 2015 
and therefore will help to improve the scientihc outcome of this mission. 
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Table 1: Stoichiometry 


Date 

YYYY-MM-DD 

Spacer 

[cm] 

T 

[K] 

-^C//4,5870 

% 

FcHijnb 

% 

FcHifibbO 

% 

CHi 

% 

2011-07-05 

n/a 

0.520 

60 

5.21 

4.26 

4.76 


2011-07-05 

n/a 

0.520 

55 

5.27 

4.26 

4.82 


2011-07-05 

n/a 

0.520 

50 

5.29 

4.28 

4.89 


2011-07-05 

n/a 

0.520 

45 

5.23 

4.61 

5.17 


2011-07-05 

n/a 

0.520 

40 

5.35 

4.40 

5.06 


2011-07-05 

n/a 

0.520 





4.86T0.41 

2011-07-06 

n/a 

0.520 

60 

8.92 

9.20 

9.45 


2011-07-06 

n/a 

0.520 

55 

9.00 

9.29 

9.63 


2011-07-06 

n/a 

0.520 

50 

9.08 

9.31 

9.64 


2011-07-06 

n/a 

0.520 

45 

9.13 

9.37 

9.84 


2011-07-06 

n/a 

0.520 





9.32±0.28 

2011-07-07 

n/a 

0.520 

60 

16.04 

16.21 

16.54 


2011-07-07 

n/a 

0.520 

55 

16.16 

16.33 

16.80 


2011-07-07 

n/a 

0.520 





16.35±0.28 

2013-06-24 

KBr 

0.025 

70 

85.08 

84.71 

88.06 


2013-06-24 

KBr 

0.025 

60 

85.89 

82.77 

87.58 


2013-06-24 

KBr 

0.025 

55 

86.43 

83.68 

86.82 


2013-06-24 

KBr 

0.025 





85.67±1.77 

2013-06-25 

KBr 

0.025 

80 

91.30 

89.10 

96.13 


2013-06-25 

KBr 

0.025 

70 

90.86 

90.47 

88.91 


2013-06-25 

KBr 

0.025 

60 

91.73 

89.15 

96.12 
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2013-06-25 

KBr 

0.025 

55 

91.74 

89.71 

96.68 


2013-06-25 

KBr 

0.025 

50 

91.31 

90.37 

97.09 


2013-06-25 

KBr 

0.025 





92.04±2.93 

2013-06-26 

KBr 

0.025 

70 

71.01 

68.00 

75.16 


2013-06-26 

KBr 

0.025 

60 

71.28 

70.05 

75.76 


2013-06-26 

KBr 

0.025 





71.88±3.01 

2013-07-10 

n/a 

0.522 

60 

20.89 

21.98 

21.84 


2013-07-10 

n/a 

0.522 





21.57±0.59 

2013-07-17 

n/a 

0.522 

60 

5.46 

5.76 

4.51 


2013-07-17 

n/a 

0.522 

55 

5.28 

5.80 

5.48 


2013-07-17 

n/a 

0.522 

50 

5.26 

5.74 

5.59 


2013-07-17 

n/a 

0.522 

45 

5.38 

5.91 

5.38 


2013-07-17 

n/a 

0.522 

42 

5.32 

5.93 

5.82 


2013-07-17 

n/a 

0.522 

40 

5.45 

6.61 

6.68 


2013-07-17 

n/a 

0.522 





5.63 ±0.49 

2013-07-18 

n/a 

0.522 

60 

2.47 

2.05 

1.80 


2013-07-18 

n/a 

0.522 

55 

2.51 

2.10 

1.55 


2013-07-18 

n/a 

0.522 

50 

2.33 

2.26 

2.47 


2013-07-18 

n/a 

0.522 

45 

2.23 

2.70 

3.60 


2013-07-18 

n/a 

0.522 

42 

2.46 

2.92 

3.06 


2013-07-18 

n/a 

0.522 

40 

2.41 

3.02 

3.56 


2013-07-18 

n/a 

0.522 

38 

2.43 

3.15 

3.63 


2013-07-18 

n/a 

0.522 

36 

2.45 

3.14 

3.17 


2013-07-18 

n/a 

0.522 





2.64 ±0.56 

2013-07-19 

n/a 

0.522 

60 

11.62 

12.44 

12.25 
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2013-07-19 

n/a 

0.522 

55 

11.63 

12.77 

13.24 


2013-07-19 

n/a 

0.522 

50 

11.79 

13.16 

13.46 


2013-07-19 

n/a 

0.522 





12.48 ±0.71 

2013-08-16 

CaFs 

0.011 

70 

85.64 

78.96 

83.40 


2013-08-16 

CaFs 

0.011 

60 

85.84 

78.15 

82.11 


2013-08-16 

CaFa 

0.011 

55 

85.25 

80.55 

84.31 


2013-08-16 

CaFa 

0.011 





82.69 ±2.91 

2013-08-23 

CaFa 

0.011 

70 

70.27 

68.48 

69.71 


2013-08-23 

CaFa 

0.011 

60 

72.36 

63.67 

76.89 


2013-08-23 

CaFa 

0.011 





70.23 ±4.36 
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Table 2: Optical constants of two different CH 4 -N 2 mixtures at two different temperatures 
in the range from 5500 to 5900 cm“^. 


Wavenumber (cm 
5500.08 
5507.35 
5514.62 
5521.89 
5529.16 
5536.43 
5543.71 
5550.98 
5558.25 
5565.52 
5572.79 
5580.06 
5587.33 
5594.60 
5601.88 
5609.15 
5616.42 
5623.69 
5630.96 
5638.23 
5645.50 
5652.78 
5660.05 
5667.32 
5674.59 
5681.86 
5689.13 
5696.40 
5703.67 
5710.95 
5718.22 
5725.49 
5732.76 
5740.03 
5747.30 
5754.57 
5761.85 
5769.12 
5776.39 
5783.66 
5790.93 
5798.20 
5805.47 
5812.74 
5820.02 
5827.29 
5834.56 
5841.83 
5849.10 
5856.37 
5863.64 
5870.92 
5878.19 
5885.46 
5892.73 


92,04% CH 4 in N 2 


T = 


1,3198817 

1.3198864 

1.3198919 
1,3198977 
1,3199049 
1,3199121 
1,3199202 
1.3199250 
1.3199178 
1,3198898 
1,3198601 

1.3198488 
1,3198508 
1,3198495 
1.3198413 
1,3198380 
1,3198366 
1,3198371 
1,3198383 
1,3198406 
1.3198447 

1.3198489 
1,3198527 
1,3198560 
1,3198593 
1,3198625 
1.3198659 
1,3198695 
1,3198729 
1,3198766 
1,3198805 
1,3198850 
1.3198906 
1,3198967 
1,3199028 
1,3199103 
1,3199188 
1,3199279 
1.3199381 
1,3199441 
1,3199365 
1,3199021 
1,3198574 
1,3198355 
1.3198333 
1,3198375 
1,3198431 
1,3198498 
1,3198564 
1,3198618 
1.3198674 
1.3198733 
1,3198789 
1,3198853 

1.3198919 


k 

0.0000040 

0.0000054 

0.0000077 

0.0000111 

0.0000158 

0.0000237 

0.0000348 

0.0000542 

0.0000812 

0.0001000 

0.0000923 

0.0000746 

0.0000636 

0.0000659 

0.0000604 

0.0000527 

0.0000456 

0.0000391 

0.0000327 

0.0000269 

0.0000224 

0.0000194 

0.0000174 

0.0000159 

0.0000146 

0.0000133 

0.0000120 

0.0000115 

0.0000109 

0.0000102 

0.0000097 

0.0000092 

0.0000092 

0.0000108 

0.0000125 

0.0000153 

0.0000206 

0.0000286 

0.0000422 

0.0000655 

0.0000984 

0.0001253 

0.0001181 

0.0000899 

0.0000656 

0.0000499 

0.0000387 

0.0000310 

0.0000264 

0.0000227 

0.0000195 

0.0000176 

0.0000160 

0.0000150 

0.0000155 


1.3198855 

1.3198900 

1.3198955 

1.3199019 

1.3199099 

1.3199179 

1.3199281 

1.3199366 

1.3199321 

1.3198977 

1.3198591 

1.3198484 

1.3198535 

1.3198546 

1.3198454 

1.3198408 

1.3198388 

1.3198392 

1.3198408 

1.3198434 

1.3198473 

1.3198518 

1.3198558 

1.3198594 

1.3198626 

1.3198659 

1.3198696 

1.3198730 

1.3198761 

1.3198797 

1.3198832 

1.3198878 

1.3198932 

1.3198994 

1.3199056 

1.3199134 

1.3199224 

1.3199327 

1.3199451 

1.3199555 

1.3199521 

1.3199118 

1.3198533 

1.3198298 

1.3198310 

1.3198379 

1.3198450 

1.3198527 

1.3198596 

1.3198649 

1.3198707 

1.3198767 

1.3198822 

1.3198884 

1.3198950 


50 K 

k 

0.0000033 

0.0000044 

0.0000059 

0.0000088 

0.0000129 

0.0000204 

0.0000308 

0.0000514 

0.0000846 

0.0001113 

0.0000989 

0.0000755 

0.0000623 

0.0000666 

0.0000625 

0.0000547 

0.0000467 

0.0000392 

0.0000327 

0.0000267 

0.0000220 

0.0000187 

0.0000167 

0.0000152 

0.0000138 

0.0000124 

0.0000115 

0.0000112 

0.0000106 

0.0000098 

0.0000092 

0.0000085 

0.0000081 

0.0000092 

0.0000106 

0.0000128 

0.0000172 

0.0000245 

0.0000375 

0.0000617 

0.0001012 

0.0001395 

0.0001296 

0.0000914 

0.0000633 

0.0000468 

0.0000361 

0.0000289 

0.0000249 

0.0000213 

0.0000184 

0.0000166 

0.0000151 

0.0000140 

0.0000142 


12.48% CH 4 in N 2 
T = 60 K T = 


1.3199970 

1.3199974 

1.3199978 

1.3199983 

1.3199989 

1.3199997 
1.3200006 
1.3200017 
1.3200030 
1.3200043 
1.3200024 

1.3199948 
1.3199916 
1.3199922 
1.3199935 
1.3199925 
1.3199908 
1.3199903 
1.3199903 
1.3199906 
1.3199910 

1.3199914 
1.3199918 
1.3199922 
1.3199925 
1.3199929 
1.3199933 
1.3199937 
1.3199942 
1.3199944 

1.3199949 
1.3199954 
1.3199960 
1.3199965 
1.3199972 
1.3199980 
1.3199987 

1.3199998 
1.3200010 
1.3200024 

1.3200042 
1.3200057 

1.3200043 
1.3199944 
1.3199866 
1.3199866 
1.3199880 
1.3199890 
1.3199899 
1.3199908 

1.3199915 
1.3199922 
1.3199928 
1.3199935 
1.3199942 


k 

0,0000004 

0,0000006 

0,0000008 

0,0000010 

0.0000013 

0.0000017 

0,0000022 

0,0000032 

0,0000048 

0,0000081 

0.0000139 

0.0000155 

0,0000109 

0,0000083 

0.0000078 

0,0000093 

0,0000081 

0.0000066 

0,0000054 

0,0000045 

0,0000038 

0,0000033 

0,0000029 

0.0000025 

0,0000022 

0,0000020 

0,0000018 

0,0000017 

0,0000015 

0.0000014 

0,0000013 

0,0000013 

0,0000012 

0,0000012 

0,0000013 

0.0000014 

0,0000016 

0,0000020 

0,0000027 

0,0000038 

0,0000058 

0.0000098 

0,0000172 

0,0000220 

0,0000155 

0,0000096 

0,0000068 

0.0000052 

0.0000043 

0,0000036 

0,0000030 

0,0000026 

0,0000023 

0.0000020 

0.0000017 


1.3199971 

1.3199975 

1.3199980 
1.3199986 
1.3199992 
1.3200001 

1.3200011 
1.3200023 
1.3200039 
1.3200057 
1.3200045 
1.3199947 
1.3199911 

1.3199921 
1.3199937 
1.3199928 
1.3199908 
1.3199902 
1.3199902 
1.3199906 
1.3199910 
1.3199914 
1.3199918 

1.3199922 
1.3199925 
1.3199930 
1.3199934 
1.3199939 
1.3199943 
1.3199947 
1.3199950 
1.3199955 
1.3199961 
1.3199966 
1.3199973 

1.3199981 
1.3199990 
1.3200001 

1.3200012 
1.3200029 
1.3200051 
1.3200072 
1.3200068 
1.3199950 
1.3199850 
1.3199856 
1.3199877 
1.3199888 
1.3199899 
1.3199908 
1.3199917 

1.3199923 
1.3199930 
1.3199936 
1.3199943 


50 K 

k 

0.0000002 

0.0000003 

0.0000005 

0.0000006 

0.0000009 

0.0000012 

0.0000018 

0.0000027 

0.0000044 

0.0000078 

0.0000149 

0.0000175 

0.0000113 

0.0000083 

0.0000078 

0.0000098 

0.0000086 

0.0000068 

0.0000055 

0.0000045 

0.0000038 

0.0000032 

0.0000028 

0.0000025 

0.0000022 

0.0000019 

0.0000017 

0.0000016 

0.0000015 

0.0000014 

0.0000013 

0.0000012 

0.0000010 

0.0000010 

0.0000011 

0.0000012 

0.0000014 

0.0000017 

0.0000023 

0.0000033 

0.0000052 

0.0000093 

0.0000179 

0.0000252 

0.0000166 

0.0000095 

0.0000065 

0.0000050 

0.0000040 

0.0000034 

0.0000028 

0.0000025 

0.0000022 

0.0000019 

0.0000017 



